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Abstract 

A QCD analysis of the world data on polarized deep inelastic scattering is presented in 
next-to-leading order, including the heavy flavor Wilson coefficient in leading order in the 
fixed flavor number scheme. New parameterizations are derived for the quark and gluon 
distributions and the value olas{M^) is determined. The impact of the variation of both the 
renormalization and factorization scales on the distributions and the value of a^ is studied. 

We obtain a^^*-'(M|) = 0.1132 „'„„„_ • The ffist moments of the polarized twist-2 
"* ^ ^^ -0.0095 

parton distribution functions are calculated with correlated errors to allow for comparisons 

with results from lattice QCD simulations. Potential higher twist contributions to the 

structure function gi(x,Q^) are determined both for proton and deuteron targets. 



1 Introduction 

The short-distance structure of the nucleon spin still consists a developing field. Nucleons as 
composite fermions obtain their spin in terms of a superposition of the spins and orbital angular 
momenta of their constituents, the quarks and gluons. It came as a surprise when the European 
Muon Collaboration (EMC) published its result [1] more than 20 years ago, which showed that 
the quarks do contribute only by a small fraction to the nucleon's spin. The obvious conclusion 
was to assume that the spin of the gluons and the orbital angular momenta of all constituents 
have to account for the missing fraction. This result initiated activities worldwide both on the 
experimental and the theoretical side in order to understand this spin puzzle and, finally, the 
spin structure of the nucleon. 

Experiments performed at CERN, SLAC, DESY and JLAB [THT2] have contributed a vast 
amount of experimental data on inclusive polarized deeply inelastic lepton-nucleon scattering 
(DIS) during the last years. The main interest in measuring the short distance structure of 
polarized nucleons has somewhat moved from determining the first moments of the twist-2 
parton distributions to the extraction of their x-dependence in the measured region and their 
scaling violations due to QCD-evolution. At large enough four-momentum transfer Q^ = — g^, 
the structure function gi{x,Q'^) mainly receives twist-2 contribution^ and is related to the 
polarized twist-2 parton distributions. 

In the present paper a QCD analysis of the polarized deep-inelastic world data is performed at 
next-to-leading order (NLO). Due to a larger set of new data which has become available recently 
the present analysis extends and updates earlier investigations [^[T5 |[T7H27] . In the QCD-fit 
we determine the fiavor singlet and non-singlet contributions of the polarized parton densities 
together with the QCD-scale Aqcd with correlated errors. The measurement of the strong 
coupling constant as(M|) from polarized deep-inelastic data does not reach the same precision 
as in the unpolarized case [281431] since the measurement is based on an asymmetry and the 
present analysis is performed in NLO. However, a consistent analysis requires the determination 
of the QCD-scale Aqcd along with the parameters of the non-perturbative input distributions. 
Also, it is interesting to see which value of as(M|) is obtained in comparison to other deep- 
inelastic analyses. At a given scale Q^ the Mellin moments of the parton distribution functions 
can be calculated under some assumptions on their extrapolation outside the measured region 
towards small and very large values of the Bjorken variable x. We also analyze, to which extend 
the present data contain higher twist contributions. 

The paper is organized as follows. In Section 2 the basic formalism is lined out. The data 
analysis is described in Section 3. In Section 4 details of the fit are given and Section 5 deals 
with the error analysis. The results of the QCD analysis are presented in Section 6. In Section 7 
we determine potential higher twist contributions and Section 8 contains the conclusions. In the 
Appendix we describe the FORTRAN-code through which the polarized parton distributions and 
structure functions can be obtained for numerical analyses. 

2 Basic Formalism 

The twist-2 contributions to the spin-dependent structure function gi{x, Q"^) are given in terms 
of a Mellin convolution of the polarized singlet AS, the gluon AG and the fiavor non-singlet 



^Twist-3 contributions are connected by target mass effects, cf. [I6] 



(NS) Aqf^ densities with the corresponding Wilson coefficient functions ACf by 
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Here x is the Bjorken variable, e, denotes the charge of the ith quark flavor in units of the 
elementary charge and Nf is the number of light flavors. The scale /ij denotes the factorization 
scale which is introduced to absorb the collinear singularities into the renormalized partonic 
distribution functions. In addition to the factorization scale there is the renormalization scale 
/ir of the strong coupling constant ^^(/i^). The gluonic Wilson coefficient ACg accounts for the 
massless as well the massive contributions due to charm quark production for W^ > {2mc + Tn]^Y 
with rric. = 1.5 GeV, 02], at ffist order, g For the implementation of the Wilson coefficients in 
Mellin space we refer to |3l]. The parton densities and the Wilson coefficient functions are 
dependent on these scales and obey corresponding renormalization group equations, while the 
structure function gi{x, Q^), as a physical observable, is independent of the choice of both scales 
fif and /i^. 

The polarized singlet and non-singlet parton densities which occur in Eq. ([1]) are expressed 
by the individual quark flavor contributions as 



r 

i=l 

Agf^(z,/i5) = Aq,{z,fi})+Aq,{z,fi})-^AJ:{z,fi}) 



(2) 
(3) 



where Ag, denotes the polarized quark distribution of the ith light flavor. We will consider 
three light quark flavors in the present analysis and treat charm quark production through the 
gluon fusion process in leading order |32] in the fixed flavor number scheme for partons. As 
well known, the following flavor non-singlet combinations contribute in the case of pure photon 
exchange considering proton and neutron targets. 
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- [(Am + Am) + {Ad + Ad)] (As + As) 
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- [{Au + An) - {Ad + Ad)] . 



(4) 
(5) 



As usual one assumes 
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D 



Aqi 
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(6) 



although later highest precision measurements could even reveal a breaking of this relation. The 
distribution 



Ad- Au 



(7) 



cannot be measured in DIS data, whereas polarized Drell-Yan data are well suited for its deter- 
mination, cf. e.g. [31]. For the measurement of the polarized strangeness distribution function 



^2nd order corrections were calculated in the asymptotic range Q^ 3> m^ in Refs. [33] 



one would wish to have sufficiently precise polarized di-muon samples available. Nonetheless 
ffist flavor separations have been attempted in performing SIDIS analyses, cf. |25j . 
The running coupling constant a^ = a;s/(47r) is obtained as the solution of 



r) 



d log(/x^ 

where, in the MS-scheme, the coefficients of the /3-function are given by 
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/3o = -CA--TFNf, 
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/9i = -Cl--CATFNf-ACFTFNf , (9) 

with the color factors Ca = 3, Tf = 1/2, and Cf = 4/3. The matching of the scale Aq^^ is 
performed at Q^ = m^, ml, with rric = 1.5 GeV and mi, = 4.5 GeV. 

In the present analysis the spin-dependent structure functions gi{x, Q"^) and gi{x, Q"^) will 
be considered referring to Nf = 3 light partonic flavors, i.e. i = u,d, s. The spin-dependent 
structure function gf{x, Q^) is represented in terms of g{{x, Q"^) and 5'"(x, Q"^) using the relation 

^1 (^, Q') = \{^- ^^^) [^1 (^' ^') + <(^' Q')] ' (10) 

where ud = 0.05 ± 0.01 is the D-state wave probability for the deuteron |35j . 

The change of the parton densities with respect to the factorization scale fi'l = Q^ is described 
by the evolution equations, which read 



dAqf^{x,Q' 
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AP(Tn^ = f APgg{x,as) 2NfAPgg{x,a,)\ . . 

The symbol ® denotes the Mellin convolution 

[A^B]{x)= / dxidx26{x - XiX2)A{xi)B{x2) . (14) 

The spin-dependent coefficient functions ACf and anomalous dimensions APij are calculated 
to next-to- leading order in the MS-scheme [361438] . which we use in the present analysis. As 
seen from Eqs. ( ITT]) and ( TT2l) . the flavor non-singlet densities Aqf^ evolve independently, while 
AE and AG are coupled in the evolution. 

In order to solve the evolution equations, a Mellin transformation of Eqs. (TTT| [T2|) and the 
polarized parton densities A/(x) is being performed by calculating its A^th Mellin moment : 

M[A/](A^) = f dx x^-^Af{x) ,N>No,NeR. (15) 

Jo 



Here Nq is chosen such that the integral flT^ converges. Under this transformation the Melhn 
convolution ® turns into an ordinary product. After the transformation has been performed 
the argument N is analytically continued to the complex plane. This also requires analytic 
continuations of harmonic sums |39j , which is outlined in Refs. |10] in detail. The fundamental 
method of solving Eqs. ( ITTl) and (fT2l) is described in the literature in detail, see e.g. Refs. [371 

mills]. 

To next-to-leading order (NLO) the solution of the flavor non-singlet and singlet evolution 
equations are given by 



Ag,f^(iV,ao),(16) 



AGW";,) ) = ll^a.Um]LiN,a„a„ni-a„Um] ( tl^^^ll'], ) • (17) 

Here Pj^jg is the leading order splitting function for the quark-quark transition. APjIjg denotes 
the NLO non-singlet '— ' splitting function, a^ = as{Q'^) and Oq = o,s{Ql), with Ql being the input 
scale. The matrices C/i and L are evolution matrices, for details see Ref. ^2]. We refrain from 
applying so-called small x resummations [13| B^. since they are very sensitive to several series 
of less singular terms [121113], which are yet unknown. Furthermore, no factorization theorem 
exists for these terms through which non-perturbative and perturbative contributions can be 
separated in a well defined way. Likewise, no other evolution equation than that governing mass 
singularities exists to deal with these terms. 

Due to the factor-structure in Eqs. (fT6!) and (TT7I) . the Gaussian error propagation of the 
input density parameters can be calculated analytically, cf. [23], for the whole Q^ region. The 
covariance matrix of the parton distributions alone is completely determined by the fit to the 
data at the input scale. 

The input distributions A^^(A^, oq), AS(A^, ao) and AG{N,ao) are evolved to the scale Q^. 
An inverse Mellin-transform to x-space is then performed by a contour integral in the complex 
plane around all singularities on the real axis for x < Xq < 1, which can be written as 

A/(x) = - / dz\m [exp(20)a;-"(")A/[c(z)]] . (18) 

^ Jo 

In practice an integral along the path c{z) = ci + p[cos(0) + zsin(0)], with ci = 1.1, p > 
and (p = (3/4) TT is performed. The upper bound on p has to be chosen in accordance with 
the numerical convergence of the integral ( 1T8l) in practice. The result Af{x) for the respective 
distribution depends on the parameters of the spin-dependent parton distributions chosen at the 
input scale Ql, which are determined by a fit to the data and to Aqcd or as^M^), respectively. 

3 Data Analysis 

The QCD analysis being performed in the following is based on the spin-dependent structure 
functions Qi' '"'{x,Q'^). These structure functions are extracted from the experimental cross 
section asymmetries for longitudinally polarized leptons scattered off longitudinally polarized 
nucleons, 

All = ^1^ . (19) 



The arrows ^(<^) denote parallel (anti-parallel) relative spin orientation of the incoming lepton 
and nucleon. The structure function ratio gi/ Fi and the longitudinal virtual-photon asymmetry 
Ai are related to A\\ by 

91 1 



Fi (1+7^) 
and 



^ + (7 - ri)A2 



(20) 



Ai = ^-VA2. (21) 



with 



9i 



The asymmetry A2 is the transverse virtual-photon asymmetry and constitutes only a small 
correction to gi. Its contribution has been treated differently by various experiments as will be 
discussed below. The other variables are given by 

_ l-{l-y)e 

^ - l + 6i?(x,Q2) ' ^^^^ 

7 = -^= , 24 

'^y (25) 



l-e{l-y) ' 

4(1 -y)- -iW 
2t/2 + 4(l-y) + 7V 



2„,2 ■ 



(26) 



Here D denotes the virtual photon depolarization factor. It determines the fraction of the incom- 
ing lepton polarization transferred to the virtual photon. The variables e, 7 and t] are kinematic 
factors, M denotes the mass of the nucleon and y = {E — E')/E is a Bjorken scaling variable 
which describes the normalized energy transfer to the virtual photon, with E the incoming energy 
and E' the energy of the scattered lepton in the target rest frame. Finally, R denotes the ratio of 
the longitudinal and transverse virtual-photon absorption cross section R{x, Q"^) = <JL/aT, which 
is experimentally well determined in the kinematic region considered in the present analysis. 

In order to obtain gi{x,Q'^) the measured ratio gi/Fi has to be multiplied by the spin- 
independent structure function Fi(x,Q^): 

g,{x,Q') = (^^] (x,Q2) X F,{x,Q') . (27) 

The structure function Fi{x,Q'^) can be calculated from the structure function F2{x,Q'^) by 

For R{x,Q'^) and F2{x,Q^) parameterizations of existing measurements are available as will be 
discussed below. 

The following data sets have been used in the present analysis: the EMC proton data pQ, the 
E142 neutron data [2\, the HERMES neutron data [3j, the E154 neutron data P|ll5j, the SMC 
proton and deuteron data [5j, the E143 proton and deuteron data [6], the HERMES re-analyzed 
proton and the new deuteron data [7], the E155 deuteron data [8], the E155 proton data [9], 
the COMPASS deuteron data [H], the JLAB neutron [10], proton and deuteron data [T2l[T3] . 



and the COMPASS proton data [H] |. The number of data points with Q^>1.0 GeV^ and 
W^ >3.24 GeV^ from the different data sets are summarized in Table 1 together with the x and 
Q^ ranges of the different experiments. In order to obtain the best possible statistical accuracy 
data on Ai, gi/Fi and gi are not averaged over the different Q^ values measured within a certain 
X bin. In total 1385 data points are used. Using Ai and gi/Fi data has, in addition to the higher 
statistics, the advantage of calculating gi for all these data sets in a unique way. Furthermore, 
gi data are sometimes only published as obtained from the average of asymmetries measured at 
different Q^ values, while for the QCD analysis it is important to maintain the Q^ dependence 
of the measured quantities. 

The SLAC parameterization i^iggo |l6] is used by most of the experiments when extract- 
ing gi. At the time of the EMC experiment this parameterization was not available yet and 
R was assumed to be Q^ independent. SMC adopted a combination of -R1990 (for x > 0.12) 
and a parameterization derived by NMC |17] (for x < 0.12). In the E155 experiment a re- 
cent SLAC parameterization for R, i?(1998) |18], was used. The changes in the data caused 
by using the different i?-parameterizations, however, are not significant and stay within the 
experimental errors Q For all Ai and gi/Fi data sets entering the present QCD analysis the 
SLAC -Ri99o(x, Q^) IIS] and the NMC F2(x, (5^)-parameterization |1S] is used to perform the 
calculation of gi. The same parameterizations were used by the E154 experiment while JLAB 
applied the i?(1998) SLAC parameterization of R. 

The magnitude of A2 has been measured by SMC [50], E154 [51], E143 [6], E155x [52] and 
JLAB [To] and was found to be small. Its contribution to gi/Fi and Ai is further suppressed by 
the kinematic factors 7 and rj and could in principle be neglected to a good approximation. On 
the other hand all the measurements have shown that the A2 contribution can be approximated 
by the Wandzura-Wilczek expression [53] , which is calculated from the spin-dependent structure 
function gi{x,Q'^) assuming that twist-2 contributions are dominant according to|j 

A2{x,Q') = j!^^f^}^ [g,ix,Q')+g2{x,Q'] 



WW 7(x,(5^) [^ dz 
Fi{x,Q') 



[-gi{z,Q'). (29) 
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The E143 experiment has exploited its measurement of A2 at 29.1 GeV and used the Wandzura- 
Wilczek expression to account for A2 for the other two lower beam energies. The measurement 
of A2 by E154 and E155x was done after having published the data on Ai and was therefore not 
available for a A2 correction of Ai. While E154 neglected A2, the E155 experiment has used the 
Wandzura-Wilczek approximation throughout its data. The JLAB measurement of A2 went into 
the extraction of g^. In HERMES the A2 contribution to g^' /Ff' data has been accounted for 
by using a parameterization for A2 obtained from a fit ^2 = CM^xj \fQ^ to all available proton 
and deuteron data, [7]. For all Ai data sets used in this analysis g\ has been calculated with 
the application of the Wandzura-Wilczek correction for A2. 

The target mass corrections to the structure function g\{x,Q'^^ are given by [T6|[59] 

X gi{i,Q^) 4MV x + i r'dii 

e (1 + 4M2a;VQ2)3/2 + Q2 ^^i + 4M2xVQ2)2 j^ ^^ 



9l {X,W ) — c (^ , A^j2r„'2 lr^2\?.l2^ r\2 Ci^ I A'\^2^2 lr\2\2 / t 3\\S\-,H ) 



■^Earlier data from Ref. [l^ are not considered. 

*The EMC proton data, where the biggest impact is expected, change by a few percent only, see Ref. |18) . 

^Note that this relation holds also in the presence of quark and target mass corrections |16[[?il[55] . for non- 
forward scattering 156' , for diffractive scattering ^571 , and the gluonic contributions to heavy flavor production I58j . 
Related integral relations for twist-3 contributions and structure functions with electro-weak couplings were 
derived in Refs. [T61I551. 



4MV 2-4MV/Q2 r^ d^^ [^ d^2 ^2 



Here M denotes the nucleoli mass and 

2x 
^^ l + (l+4M2xVQ2)i/2 (2^) 

is the Nachtmann variable. We refrain to use approximations to f l5U]) in terms of power series 
in M^/Q^ since the convergence of the corresponding series is problematic as being lined out in 
Ref. [16]. 

The data sets contain both statistical and systematic errors. It is known that the systematic 
errors are partly correlated, which would lead to an overestimation of the errors when added 
in quadrature with the statistical ones, and hence to a reduction of the x^ value in the fitting 
procedure. To treat all data sets on the same footing only statistical errors were used. However, 
a relative normalization shift, A/i, between the different data sets was allowed within the normal- 
ization uncertainties, AA/",, quoted by the experiments. These normalization shifts were fitted 
once and then fixed, see Table 1. Thereby the main systematic uncertainties coming from the 
measurements of the luminosity and the beam and target polarization were taken into account. 
The normalization shift for each data set enters as an additional term in the x^^^xpression for 
the fit which then reads 



X 



j=i 



(AM)2 



E ^^^9u!h - 9\T)iCr')"'i^^9u!k ' 9\T) 



h,k=l 



where the sums run over all data sets and in each data set over all data points. The covari- 
ance matrices Cj are diagonal for each experiment except for the case of the HERMES data [7] . 
The statistical errors, and consequently the covariance matrices Cj, have been rescaled by the 
normalization factors A/i. The minimization of the x^ value above to determine the best param- 
eterization of the polarized parton distributions is performed using the program MINUIT [50] . 
Only fits giving a positive definite covariance matrix at the end have been accepted in order to 
be able to calculate the fully correlated la statistical error bands. 

4 Parameterization of the Polarized Parton Distributions 

The shape chosen for the parameterization of the polarized parton distributions A/j(x, Q^) 
in x-space at the input scale Ql is : 



xAfiix, Ql) = ViA.x'^^il - xf' [1 + pix'^ + 7ixj . (32) 

The normalization constant Ai, being given by 

A"' = (1 + 7. ^ ^""l ^ J B{a,,b, + l)+p^B U + l,b, + l\ , (33) 



is calculated such that 



-i 
r]i= I dxAfi{x,Ql) (34) 







is the first moment of A/j(x,(5o) ^^ ^^^ input scale. Here, B{a,b) is tlie Euler Beta-function 
being related to the F-function by B{a, h) = T{a)T{b)/T{a + h). The choice of the shape (!32|) is 
apphed in various QCD analyses of unpolarized data, see e.g. Ref. [61]. The term x"' controls 
the behavior of the parton density at low and (1 — x)''* that at large values of x, respectively. 
The remaining polynomial factor accounts for potential degrees of freedom at medium x. At 
the same time not too many parameters can be sufficiently constrained by the data available at 
present. In particular the parameters pj, which play a role in unpolarized analyses, [28], are set 
to zero, cf. also [23|[25]. 

If the QCD evolution equations are solved in Mellin space as described in Section [2] a Mellin 
transformation of the polarized parton density A/(x, Q^) is performed for complex arguments 
N : 



M[A/,(x, g2)](iv) = / x^-it/xA/,(x,Qg 

^0 



= r/,A, fl + 7,^^^^)s(iV-l + a„6, + l) . (35) 

V N + ai + hiJ 

In the present analysis we assume (approximate) SU(3) flavor symmetry the sea quark distribu- 
tion is given by 

Ag,(x, g2) = Am(x, Q2) = a5(x, Q^) = As{x, Q^) = As{x, Q^) , (36) 

AE(x,Q2) = AM,(x,g2) + Ad,(x,Q') + 6Ag,(x,Q2) . (37) 

We refer to the inclusive polarized DIS World Data only. A breaking of the flavor symmetry 
also for the light (sea) quarks in the polarized case is probable and has been clearly observed 
in the unpolarized case, cf. e.g. [31]. In the polarized case first attempts have been made to 
determine the individual sea quark distributions, cf. Refs. [25l|62]. We consider the evolution of 
the complete light polarized sea. 

Under the above assumption four spin-dependent parton densities have to be determined in 
the QCD analysis. They are chosen to be: Am^(x,(5^), A(i„(x,(5^), Ag^(x,(5^) and AG(x, Q^). 
As seen from Eq. (!32|) . each spin-dependent density contains five parameters which gives a total 
of 20 for all four. It has been found that, in order to meet the quality of the available data and 
the reliability of the fitting procedure, this large number of free parameters has to be reduced, 
which is discussed below. 

The first moments of the polarized valence distributions Am^ and A(i„, rju^ and 77^^, can be 
fixed exploiting the knowledge of the parameters F and D as measured in neutron and hyperon 
;5-decays according to the relations: 

riu^-r]d. = F + D, (38) 

Vu^+Vd. = SF-D. (39) 

A re-evaluation of F and D was performed on the basis of updated /3-decay constants [63] 
leading to 

F = 0.464 ± 0.008 and D = 0.806 ± 0.008 , (40) 

and consequently to 

riu^ = 0.928 ± 0.014 and rj^^ = -0.342 ± 0.018 . (41) 

In order to compensate for the present insufficient accuracy of the data, a certain number of 
parameters is set to zero after this was thoroughly suggested by initial fits. This applies to 



Pu„ = Pdv = 0, 7g^ = /?g^ = 0, and 7g = Pg = 0. The number of parameters to be fitted for each 
polarized parton density is reduced to three, i.e. to 12 in total. In addition the QCD scale Aqcd 
is fitted. 

In the analysis it turns out that the four parameters 7^^, 7^^, 6,^, and he have very large 
uncertainties. The precision of the data is not high enough to constrain these parameters suf- 
ficiently. Altering them within these uncertainties does not lead to a significant change of 
X^. These four parameters were therefore fixed. The first two of them were fixed at their 
values obtained in the initial fitting pass, 7^^ = 27.64 and 7^^ = 44.26. In fixing the high- 
X slopes ho and hq^ a relation was adopted as derived from the unpolarized parton densities, 
^g^^clpol) = hqJhGiyiYqioX) = 1.44. Fitting with this constraint led to the following choice: 
be = 5.61 and bg^ = 8.08, see e.g. Ref. |6l]. 

A second relation was adopted to constrain the low-x behavior of the spin-dependent gluon 
density with respect to the low-x behavior of the spin-dependent sea-quark distribution by 
O'G = (^q^ ~^ C with C = 1. This relation, together with the relation for the high-x slopes, are 
suited to establish positivity for AG and Aq^. No explicit positivity constraint has been enforced 
for Au^ and Ady. 

5 Determination of the Errors 



5.1 Calculation of Statistical Errors 

The evolved polarized parton densities and structure functions are functions of the input densi- 
ties. Let Af{x,Q'^]Pi\i=i^k) be the evolved polarized density at the scale Q^ depending on the 
parameters Pi\i=i^k- Then its correlated statistical error as given by Gaussian error propagation 
is 

s^ fdAfdAf 



where cov{pi,pj) are the elements of the covariance matrix determined in the QCD analysis. 
The gradients dAf /dpi at the input scale Ql can be calculated analytically. Their values at Q^ 
are evaluated through the evolution and can then be used to compute the errors according to 
Eq. ( H2|) . As shown in Section [21 the covariance matrix is completely determined by the fit at 
the input scale and does not change when the evolution is done in Mellin space. That means it 
can be used at any scale Q^. For the expressions of the gradients at the input scale see Ref. |23j . 
Apart from statistical errors the data are also subject to systematic uncertainties which are 
even partly correlated. These correlations are not published by the experiments and can hence 
not be taken into account here. In the following the influence of experimental and theoretical 
systematic uncertainties will be investigated. 

5.2 Determination of Experimental Systematic Uncertainties 

The experimental systematic uncertainties were estimated from the following sources : 

1. The variation of the data within their experimental systematic uncertainties. 

The procedure used to obtain the contribution from the experimental systematic uncer- 
tainties consists in shifting each data set by icr^ysi while leaving the other data sets at their 
central values and looking at how much the polarized distributions change. The extreme 
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changes from the 'central' distribution were taken as the systematic uncertainties. This 
was done for each of the 19 data sets used separately and, finally, the 19 contributions were 
added in quadrature to obtain the total contribution. 

2. The variation of the data within the upper and lower limits of the NMC F2 parameterization [39]. 
When calculating gi{x,Q'^) from the asymmetry data, see Section [3l the NMC F2- 
parameterization was used at its upper and lower limit to determine the changes in the 
polarized distributions compared to the 'central' curve. The extreme deviations from that 
curve were taken as the systematic uncertainties arising from the F2 parameterization. 



3. The variation of the data within the uncertainty of the R parameterization [IB]. 

When calculating gi{x,Q'^) from the asjTiimetry data, see Section [3l the SLAC -R1990 pa- 
rameterization was used at its uncertainty limits, see Ref. [16], to determine the changes 
in the spin-dependent distributions compared to the 'central' curve. The maximal de- 
viations from that curve were taken as the systematic uncertainties arising from the R 
parameterization. 

Finally, the contributions from all three sources were added in quadrature to obtain the total 
contribution at each value of x, which is shown as hatched error bands in the Figures below. 

5.3 Determination of Theoretical Systematic Uncertainties 

The theoretical systematic uncertainties were estimated from the following sources : 

1. The variation of the factorization and the renormalization scale by a factor of 2. 

2. An additional variation of AL^^. We varied Aq^^ by ±30 MeV, which corresponds to a 
variation of as{M'^) by ±0.002, being a typical error in individual precision measurements, 
in addition to the error of as(M|) being determined in the present fit. 

3. The variation of 7]^^ and 77^^ within the errors of the parameters F and D, see Section 4, 
while keeping qa/ Qv = F + D constant. 

4. The variation of the parameterizations at the input scale Q^. 

Two cases are considered: first, the values of the parameter p = —1.0 for Au^ and Arf^, 
are added in the fit. Second, the values of the parameter 7 = 17.64 for Au^ and 24.26 for 
Adv, respectively, for the same densities are changed compared to the values used for the 
central curve. 

5. The variation of the standard input scale: from Ql = 4.0 GeV^ to 2.0 GeV^ and 6.0 GeV^. 

For all items always the extreme deviations from the 'central' curve were taken as the contri- 
bution to the theoretical systematic uncertainty from that source. Finally, all different contribu- 
tions were added in quadrature to get the total contribution at each value of x which is shown 
as the second hatched error bands in the Figures below. 
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6 Results of the QCD Analysis 

6.1 The parton distributions 

In the NLO QCD fit we covered the polarized world deep-inelastic data with Q^ > 1 GeV^, W'^ > 
3.24 GeV^, with a value of x^ /NDF = 1.12. We have chosen a value UP = 9.3 corresponding 
to the 7+1 parameter fit in the MINUIT-procedure. The parton distributions are parameterized 
at Qq = 4 GeV^. In Table 2 the values of the fit parameters are summarized. The covariance 
matrix of the fit is given in Table 3. Gaussian error propagation allows to derive the error bands 
due to the parton densities for the various polarized observables. 

In Figure 1 the polarized momentum distributions xAfi{x) are presented at the scale Q^. We 
compare with other analyses [T8 | [T9 | |2 H I25] . The xAu^ distribution is slightly lowered if compared 
to our previous analysis p3]. The most important change concerns the gluon distribution xAG, 
which is lowered by about a factor of two relative to the results of Ref. [23]. Comparing to 
the results of other analyses xAu^ turns out to be lower in a wider range. In case of the 
xAii^-distribution all fits widely agree within the la error band. This also applies for the xAq- 
distribution. For the polarized gluon distribution the agreement of the different fits at larger 
values of x agree within the la error, while below x ~ 0.02 the fits p ^ [T9 | [2T| yield slightly 
higher values. We indicated the positivity constraints referring to the unpolarized PDFs [30] , 
which are obeyed. 

In Figures 2 and 3 we analyze the systematic errors in more detail, cf. Sections 5.2 and 5.3. 
For the gluon distribution function xAG{x, Qq), Figure 2, the experimental systematic errors are 
lower than the statistical errors, but are still significant. The combined theoretical systematic 
effects at NLO amount to larger values than the experimental ones. About half of the error is 
due to the uncertainty in Aqcd- Clearly, in future analyses based on NNLO QCD evolution this 
error will diminish. The effect of the experimental and theoretical systematic errors in case of 
the singlet distribution xAE(x, Qo) is similar. The errors are smaller if compared to the gluon 
distribution. The singlet distribution a.t Qq = 4 GeV^ is negative (within the la errors) for 
X < 2 ■ 10~^ and turns to positive values for x > 4 ■ 10~^. All the fits p ^ [T9 | l2T] lie inside the la 
error band. In the medium x-range the DSSV distribution [2S] yields somewhat larger values. 

In Figure 4 we compare the fit results for the structure functions g^{x,Q'^),gf{x,Q^) and 
gi{x, Q"^) at Q^ = 5 GeV^ with the data, cf. Table 1, to illustrate the fit quality for the different 
targets as an example. Furthermore, also the results of the GRSV [19] and AAC p^ analyses 
are shown. Overall a good agreement is obtained. 

A further illustration of the fit quality is presented in Figure 5. Here the data for (^f (x, Q^) 
are compared to the fit including, the statistical errors. We also show the fit results of AAC, 
GRSV, and LSS [T8l[T9l|2T]. Within the error bands the data agree well with the fit. In the 
lowest x-bins the fiuctuation is somewhat larger. In some cases the EMC-data [1] lay outside 
the la error range. 

6.2 Aqcd and a^Mf) 

The NLO QCD-analysis of the polarized world deep-inelastic data requires the fit of Aqcd along 
with the parameters of the non-perturbative parton distributions at the scale Qq. As outlined 
in Table 3 the QCD-scale is correlated to all other parameters and in particular also to rjc- Due 
to this analyses in which as(M|) or Aqcd is imported from 3rd sources may suffer significant 
biases, i.e. a too large value of as(M|)(AQCD) leads to a too small gluon distribution, aside of 
other effects. Despite of various precision measurements of the strong coupling constant based 
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on theoretical NNLO (and partially even higher) precision, a thorough agreement on the value of 
as^M"^) has not yet been reached, cf. e.g. [65]. Due to this Aqcd is determined in this analysis. 
We refer to Aq(-.j-)' ^~ as the NLO value for 4 active flavors and obtain 

Ag),j3 = 243 ± 62(exp) MeV . (43) 

In an earlier analysis [23] the values 

AgJ,j^ = 235 ± 53(exp) MeV ISET = 3 (44) 

Ag),j^ = 240 ± 60(exp) MeV ISET = 4 (45) 

were found, for comparison, slightly depending on some assumptions in the fit. The variation of 
the factorization and renormalization scales filthy a. factor of 1/2 and 2, respectively, yields 

Agio = 243 ±62 (exp) '^ (FS) +g^ (RS) MeV. (46) 

Here we excluded values fi'j,^ < 1 GeV^, unlike in Ref. [,23j, since at scales lower than 1 GeV^ 
the perturbative description cannot be considered reliable anymore. 
Correspondingly, for q;s(M|) one obtains 

/.^2N niioo +0.0043 , , -0.0029 ,^„, +0.0032 ,^^, , ,^. 

a,(M|) = 0.1132 _Q_QQ^^ (exp) ^^ ^^^^ (FS) _^^^^^ (RS) , (47) 

with combined errors of 

a.(M|) = 0.1132 +°;°°^^ . (48) 

Due to the NLO analysis the factorization- and renormalization scale uncertainties are still 
dominant. The values are well compatible with recent determinations of the strong coupling 
constant at NNLO and N^LO from deep-inelastic data : 

«,(M|) = 0.1134 +|];J]U2? NNLO [28] (49) 

«,(M|) = 0.1141 +2;j;[|^0 N^LO [28] (50) 

a^(M|) = 0.1135 + 0.0014 NNLO, FFS [31] (51) 

a,(M|) = 0.1129 + 0.0014 NNLO, BSMN [31] (52) 

a,(M|) = 0.1124 + 0.0020 NNLO, dyn. approach [29] (53) 

as{Ml) = 0.1158 + 0.0035 NNLO, stand, approach [29] (54) 

a,(M|) = 0.1171 + 0.0014 NNLO [30] (55) 



More recent unpolarized NNLO fits, including the combined HERA data [66], yield 

a,(M|) = 0.1147 + 0.0012 NNLO [67] (56) 

a,(M|) = 0.1145 + 0.0042 NNLO, (preliminary) [68] (57) 

The central value of the present fit ( HHI) does well compare with the above values. They are 
located below the present weighted average of as(M|) measurements [65] of 

a,(M|) = 0.1184 + 0.0007. (58) 
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The error given in (1581) does not include any yet unknown relative systematics between the 
different classes of the same type of measurement. 

We would like to mention that recent determinations of as(M|) using event shape moments 
for high energy e^e~ annihilation data from PETRA and LEP including power corrections the 
following values were obtained : 

a^(M|) = 0.1135 ± 0.0002 (exp)± 0.005 (fii)± 0.0009 (pert) NNLO [69] (59) 
as{Ml) = 0.1153 ± 0.0017 (exp)± 0.0023 (th) NNLO [70] (60) 

Also these measurements of as(M|) yield low values. They show that the results obtained ana- 
lyzing deep-inelastic data do not form a special case. The systematics of the different extractions 
of as(M|) has to be understood in more detail in the future. 

Fit results from previous polarized analyses like [T5l[l71[26] were summarized in [23]. In 
Figure 6 we compare recent determinations at NNLO and N^LO for unpolarized and at NLO 
for polarized deep-inelastic scattering. 

6.3 Moments of Polarized Parton Distributions 

We calculate the lowest moments of the polarized parton densities 

{f{x)U = [ dxx'^Afix) , (61) 

Jo 

where Af{x) denote the different polarized (number) density distributions. The moments n = 
0, ...,3 are given in Table 4. The behaviour of these distributions outside the kinematic range 
in which the fit is performed bear uncertainties, which are difficult to predict for these non- 
perturbative quantities. U Instead of presenting necessarily uncertain models for this range, we 
compute the respective part of the moments for values x < 0.005 and x > 0.75 extrapolating 
the present distributions to the range x e [0, 1]. 

The zeroth moments of the polarized quark- and gluon distributions as well as the contribu- 
tions due to the quark- and gluon angular momenta, Lg and Lg, constitute the nucleon spin 



1 _ 1 

2 ~ 2 



(AS(x))o + {AG{x))o + Lg + L,. (62) 



We obtain 



{AJ:{x,QI))o = 0.193 ±0.075 (63) 

{AG{x,Ql))o = 0.462 ±0.430, (64) 



at Ql = 4GeV2. For Eq. ([62]) this yields 

1 



2 



(0.555 ± 0.436) + Lg + Lg . (65) 



The error of AG is clearly dominant. 

The results given in Table 4 can be compared to ab initio calculations of these moments 
in Lattice Gauge Theory. There we also compare with the values obtained in our previous 



^We remind the failure in predicting the lower x behaviour of the structure function F2{x,Q^) prior the 
HERA measurements until 1992, which assumed a slightly falling or constant behaviour below x ~ 10^^, whereas 
a strong rise was measured at HERA. 
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analysis [23] . In the present analysis the first moments {n = 0) of the polarized valence quark 
distributions are determined by the values of F and D and are fixed in the fit. Comparing to 
the results of [23], ISET = 31} slightly larger values are obtained for the moments of Am„(x, Q"^) 
and slightly lower values for Adv{x, Q^) and Ag(x, Q'^). A very significant change is obtained for 
the moments of the polarized gluon density, where the moments reduced by a factor of about 
two comparing to [23]. Although being positive, the latter moments are now compatible with 
zero in the la errors. 

First lattice results for the moments n = 0, 1, 2 of the polarized quark distributions were given 
about ten years ago. Meanwhile many systematic effects in the simulation have been improved 
further. Still there are differences in the different simulations. Rather aiming on a detailed 
comparison with the values in Table 4 we give a brief summary of the current status. A recent 
survey has been presented in [71]. Lattice results on (Au — A(i)o for m^ = 0.029...0.48GeV^ 
were given in [72^E5] by the BGR, RBC, LHPC, ETMC, QCDSF-collaborations using dynamical 
quarks. Most of the values are yet below the experimental value. The QCDSF collaboration [75] 
performed simulations at m^ = 170 MeV and obtained 

{Au{x) - Ad{x))o = 1.17 ±0.05, (66) 

(Au(x) - Ad(x))o , , 

)— 4^ ^ , ° = 0.47 ±0.02. 67 

For the first moment the following values were determined 



(An-Arf)i =0.271 ±0.040, 


m^ = 493 MeV, 


[72; 


(68) 


= 0.252 ±0.020, 


m^ = 352 MeV, 


74; 


(69) 


which are larger than the value 








(AM-Arf)i = 


0.190 ±0.008 




(70) 



determined in the present analysis. 

Results on the second moment were given in [71] 



to be compared to 



{Au-Ad)2 = 0.083 ±0.012, m^ = 352MeV, (71) 

(72) 



{Au-Ad)2 = 0.063 ±0.004 (73) 



obtained in this analysis. 

Results of older lattice simulations [76] were discussed in [23] previously. The values given 
above for the 1st and 2nd moment are based on theoretically much improved simulations, if com- 
pared to early investigations [76]. Yet, the moments obtained yield similar values. Comparing 
the lattice results with the results obtained in QCD-fits to the polarized deep-inelastic world 
data one observes a similar trend of values but not yet agreement. 



^Thc NLO results for ISET = 4 are quite similar. 
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7 Higher Twist 

So far we have applied the twist-2 approximation at NLO to describe the spin-dependent struc- 
ture function gi. As the data may contain contributions from higher twist (HT) it is of interest 
to look for possible effects of such contributions. A thorough description of higher twist anoma- 
lous dimensions and Wilson coefficients to NLO is still missing, even for the twist-4 contribu- 
tions. Therefore we will perform a purely phenomenological analysis. Similar to the approach of 
Ref. [77] for the structure function F2, where a higher twist term is parameterized by the ansatz 

/^W(^,Q') = ^, (74) 

and used multiplicative to the leading twist (LT) contribution, f7i(x,Q^) is described by 

9{{x, Q') = gl^A^, Q') [I + h'HT,mi^, Q')] , (75) 

where i = p,d, n. Likewise onw may describe an additive higher twist term 

glix,Q') = glLAx,Q') + h^nrA^^Q') ■ (76) 

This approach has to be handled with great care, since the coefficients Ci{x) are actually also 
Q^ dependent. They consist of a combination of various terms which exhibit different scaling 
violations. The relation to Aqcd is completely masked here. Moreover, higher twist contributions 
should have a flavor- dependence and are not expected to be the same in case of polarized and 
unpolarized deep-inelastic scattering. 

The kinematic x-range being covered by experiment is divided into 5 bins and the coefficient 
Ci{x) has to be determined in each bin and for each target. The resulting coefficients for the 
proton and the deuteron target, Cp{x) and Cd{x), are summarized in Table 5 both for the 
multiplicative and additive cases. Due to the enlarged number of parameters being fitted the 
value of UP = 11.5 is used. In the additive and multiplicative cases we observe quite comparable 
patterns. We prefer the additive case, since the twist-2 scaling violations of gi{x,Q'^) do not 
influence Cp^d,n{^)- Here, in most of the bins the result is 1 a compatible with zero, except 
of one case for Cp{x). CdA^) ^^^ ^^^ 1 ^ compatible with zero. The behaviour in case of the 
deuteron and neutron is rather flat, while for the proton a slight structure, yet with large errors, 
is indicated, cf. Figure [71 In the flt determining also the higher twist parameters the value 
X^/NDF for the CLAS data amounts to 1.12. The present data are not yet precise enough 
to undoubtedly reveal higher twist contributions in the range Q^ > 1 GeV^ and a NLO QCD 
analysis can be carried out in the leading twist approximation. We do not conflrm the results of 
Ref. [22]. Unlike the case for the large x valence quark region, in which dynamical higher twist 
terms are extracted consistently in the unpolarized case, cf. Refs. [2Hll7Sl[7n], the situation is 
more involved for the lower x-region. The dynamics is clearly different in both these domains 
due to the contributing parton species. As has been shown in [80], different power corrections 
cancel each other in the small x region. 

8 Conclusions 

A QCD analysis of the polarized deep-inelastic world data has been performed at NLO, including 
the effects of charm production to flrst order. We derived a parameterization for the polarized 
parton distributions and Aq^^ with the error correlations between the fltted parameters applying 

15 



the x^-inethod. Detailed comparisons have been performed with recent parameterizations 
[T9l[2Tl[25]. The present data are not accurate enough to determine all the shape parameters at 
a sufficient accuracy. Due to this some of the parameters have to be fixed after an initial phase 
of the analysis to a model- value. If compared to our previous analysis [23] the more recent data 
lead to a smaller gluon distribution, which is for a wide region of x compatible with zero within 
the la error. We determined both the experimental and theoretical systematic effects. Both 
the central values of the parton densities and their la error are made available in form of a 
numerical parameterization in the range x G [10~^, 1], Q^ G [1,10^] GeV^. These distributions 
can be used for polarized hard-scattering processes at hadron- and lepton-nucleon colliders for 
various observables, including error propagation w.r.t. the accuracy of the parton densities. The 
implementation in terms of grid-interpolation is well suited also for Monte Carlo simulations. 

The QCD-scale was determined by Aq^J,^ = 243 ± 62 (exp) _„„ (th) MeV, corresponding 
to a^^'^{M'^) = 0.1132 _r, nnqr . The central value is well compatible with other measurements, 



cf . [28l[29l[3Tl[67] - [7()] . The errors are still rather large, also because of the scale variation uncer- 
tainties at NLO. Nonetheless the correlated determination of as(M|) with the parton densities 
is of importance to avoid biases in particular w.r.t. to the size of the gluon distribution function. 
We also determined potential higher twist contributions, which were found to be compatible 
with zero in the whole kinematic range within the present errors for C„,d(x) and also for Cp{x), 
except for one of the bins. Based on the results of the present analysis we computed the lowest 
moments of the individual twist-2 parton densities. For the lowest moment (l/2)(AE(a;))o + 
{AG{x))o we obtain 0.555 ± 0.436 at Qq = 4 GeV^, which is well compatible with the nucleon 
spin 1/2 even without angular momentum contributions. However, the error is dominated by 
that of the polarized gluon distribution. The moments derived may be compared to upcoming 
lattice simulations based on dynamical quarks of the corresponding operator matrix elements. 
The present results are not yet in agreement, although the tendency of values is visible. Runs 
at even smaller values of 171.,^ seem to be necessary. 
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9 Appendix: The FORTRAN-code for the parton densities 
and their errors 

A fast FORTRAN program is available to represent the polarized parton densities xAm^(x, Q^), 
x/S.d^{x, Q'^), xAG{x, Q'^), and xAq{x, Q"^), as well as the polarized structure functions xg^{x, Q"^) 
and xg^{x, Q"^) at NLO in the MS-scheme together with the parameterizations of their la errors. 
The following ranges in x and Q^ are covered: 



10-9 <x<l 



1 GeV2 < g2 < 106 GeV2. 



The polarized distributions are the result of a fit to the world data on spin asymmetries, i.e. A^'""' 
or gi/Ff'^' , as described above. The SUBROUTINE POLPDF returns the values of the polarized 
distributions, always multiplied by x, at a given point in x and Q^ by interpolating the data on 
specified grids. The interpolation in x is done by cubic splines and in Q^ by a linear interpolation 

The parton distributions are evaluated by 

SUBROUTINE POLPDF (ISET, X, Q2, UV, DUV, DV, DDV, GL, DGL, SEA, DSEA, 

G1P,DG1P,G1N,DG1N), 

with ISET = 1. All non-integer variables are of the type REAL*8. The calling routine has to 
contain the COMMON/INTINI/ IINI. Before the first call to SUBROUTINE POLPDF the initialization 
is set by IINI = 0. 

The parameters X, Q^ [GeV^] are x and Q^. The momentum densities of the polarized up- 
and down valence quarks, gluons and the sea quarks are UV , DV , GL , SEA, with SEA = xAug = 
xAdg = xAu = xAd = xAs = xAs. Correspondingly, DUV is the Icr error of UV etc. and GIP 
and GIN are the values of the electromagnetic structure functions xgf and xg"^. 

The programme example. f reads the data-grid qcd_nlo_905_0.grid and is compiled using 
gf ortran at a LINUX-system. The test-code produces the test-output for the structure-functions 
xglp, xgln, xgld and their lo" errors dxglp, dxgln, dxgld : 



* X, 


,q2,glp,dglp,gln,dgln,gld,dg 


Id 




















0, 


, 100000 


4.000000 


0, 


.027274 


0, 


.001453 


-0, 


.011366 





.001468 


0, 


.007358 


0, 


.000955 


0, 


.200000 


4.000000 


0, 


.043553 


0, 


.001347 


-0, 


.007264 





.001895 


0, 


.016784 


0, 


.001075 


0, 


.300000 


4.000000 


0, 


.052548 


0, 


.001146 


-0, 


.004617 





.001947 


0, 


.022168 


0, 


.001045 


0, 


.400000 


4.000000 


0, 


.054270 


0, 


.001137 





.000329 





.002212 


0, 


.025252 


0, 


.001150 


0, 


.500000 


4.000000 


0, 


.046195 


0, 


.001369 





.003331 





.002272 


0, 


.022906 


0, 


.001227 


0, 


.600000 


4.000000 


0, 


.034999 


0, 


.001821 





.003411 





.001920 


0, 


.017764 


0, 


.001224 


0, 


.700000 


4.000000 


0, 


.022119 


0, 


.001858 





.002789 





.001254 


0, 


.011520 


0, 


.001037 


0, 


.800000 


4.000000 


0, 


.011919 


0, 


.001254 





.001508 





.000557 


0, 


.006210 


0, 


.000634 


0, 


.900000 


4.000000 


0, 


.006516 


0, 


.000383 





.000476 





.000115 


0, 


. 003234 


0, 


.000185 


0, 


.950000 


4.000000 


0, 


.005713 


0, 


.000087 





.000252 





.000023 


0, 


.002759 


0, 


.000042 



The program can be received on request via e-mail to Johannes.Bluemlein@desy.de or 
Helmut.Boettcher@desy.de or from http : //www-zeuthen . desy . de/~blumlein , 



^We thank S. Kumano and M. Miyama of the AAC-coUaboration for aUowing us to use their interpolation 
routines. 



17 



10 Tables 



Experiment 


x-range 


Q^-range 

[GeV^] 


data points 


M 


type 


# 


E143(p) [6j 
HERMES (p) [7] 
E155(p) [9] 
SMC(p) [5] 
EMC(p) [1] 
CLASl(p) I2| 
CLAS2(p) 13] 
COMPASS(p) [14] 


0.027 - 0.749 
0.026 - 0.731 
0.015 - 0.750 
0.004 - 0.484 
0.015 - 0.466 
0.125-0.575 
0.292 - 0.592 
0.005 - 0.568 


1.17-9.52 
1.12- 14.29 
1.22 - 34.72 
1.14- 72.10 
3.50 - 29.5 
1.10-4.16 
1.01-4.96 
1.10-62.10 


9i/Fi 
Ax 
Ax 
Ax 

9i/Fi 
Ax 


82 
37 
24 
59 
10 
10 
191 
15 


0.963 
0.970 
1.003 
0.960 
0.964 
1.010 
1.030 
0.955 


proton 








428 




E143(d) [6] 
HERMES (d) [1] 
E155(d) [8] 
SMC(d) [5] 
COMPASS (d) 11] 
CLASl(d) 12] 
CLAS2(d) PJ 


0.027 - 0.749 
0.026 - 0.731 
0.015 - 0.750 
0.004 - 0.483 
0.005 - 0.566 
0.125-0.575 
0.298 - 0.636 


1.17-9.52 
1.12- 14.29 
1.22 - 34.79 
1.14-71.76 
1.10-55.30 
1.01-4.16 
1.01-4.16 


9i/Fx 
Ax 

9i/Fx 
Ax 
Ax 
Ax 

9i/Fx 


82 
37 
24 
65 
15 
10 
662 


0.960 
0.970 
0.979 
0.998 
0.952 
1.003 
1.014 


deuteron 








895 




E142(n) [2] 
HERMES (n) [3] 
E154(n) 4]/ [15 
JLAB(n) pii 


0.035 - 0.466 
0.033 - 0.464 
0.017-0.564 
0.330 - 0.600 


1.10-5.50 

1.22 - 5.25 

1.20 - 15.00 

2.71-4.83 


Ax 
9i 
91 
91 


33 

9 

17 
3 


0.989 
0.970 
0.980 
1.000 


neutron 








62 




total 








1385 





Table 1: Number of data points on Ax, gx/Fx or gx for Q^ > 1.0 GeV^ and W'^ > 3.24 GeV^ used in 
the present QCD analysis. For each experiment are given the x and Q^ ranges, the type of quantity 
measured, the number of data points for each given target, and the fitted normalization shifts Mi 
(see text). 



Au^ 


V 


0.928 (fixed) 


Aq, 


V 


-0.066 ± 0.013 




a 


0.239 ± 0.027 




a 


0.365 ± 0.164 




b 


3.031 ± 0.178 




b 


8.080 (fixed) 




1 


27.64 (fixed) 




7 


0.0 (fixed) 


Ad^ 


V 


-0.342 (fixed) 


AG 


V 


0.462 ± 0.430 




a 


0.128 ± 0.068 




a 


OAg. + 1 




b 


4.055 ± 0.879 




b 


5.610 (fixed) 




7 


44.26 (fixed) 




7 


0.0 (fixed) 


A^^l^ = 243 ± 62 MeV 


X^/NDF = 1537/1377 = 1.12 



Table 2: Final parameter values and their statistical errors at the input scale Ql = 4.0 GeV^. 





a(4) 


'^ti„ 


bu^ 


ad. 


bd^ 


"Hsea 


^sea 


VG 


A (4) 


3.84E-3 
















a«„ 


-4.08E-4 


7.56E-4 














bu. 


-1.14E-3 


4.30E-3 


3.18E-2 












O'd^ 


2.74E-3 


-9.39E-4 


-4.43E-3 


4.60E-3 










bd^ 


2.38E-2 


-8.33E-3 


-1.03E-2 


4.50E-2 


7.72E-1 








1]sea 


1.70E-3 


-6.84E-4 


-3.60E-3 


2.26E-3 


2.11E-2 


5.62E-3 






O'sea 


-5.64E-3 


3.04E-3 


1.65E-2 


-8.25E-3 


-7.37E-2 


7.13E-4 


2.70E-2 




VG 


-1.96E-2 


8.32E-3 


4.25E-2 


-2.16E-2 


-1.67E-1 


-2.09E-2 


4.17E-2 


1.85E-1 



Table 3: The covariance matrix for the 7-1-1 parameter NLO fit based on the world asymmetry data. 
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Distribution 




Fit Results 


[23], set 3 


n 


value 


value out 








of range 




Am„ 





0.928 ±0.000 


0.158|3.3E-3 


0.926 ±0.071 




1 


0.153 ±0.004 


1.6E-4|2.7E-3 


0.163 ±0.014 




2 


0.052 ±0.002 


2.1E-3 


0.055 ± 0.006 




3 


0.023 ±0.001 


1.7E-3 


0.024 ±0.003 


Ad, 





-0.342 ± 0.000 


-0.1101-2. lE-4 


-0.341 ±0.123 




1 


-0.037 ±0.007 


-7.0E-5I-1.7E-4 


-0047 ±0.021 




2 


-0.010 ±0.002 


0| -1.3E-4 


-0.015 ±0.009 




3 


-0.004 ±0.001 


Ol-l.lE-4 


-0.006 ±0.005 


Au-Ad 





1.270 ±0.000 


0.267|3.5E-3 


1.267 ±0.142 




1 


0.190 ±0.008 


2.3E-4|2.8E-3 


0.210 ±0.025 




2 


0.063 ± 0.004 


2.3E-3 


0.070 ±0.011 




3 


0.027 ±0.002 


1.8E-3 


0.030 ±0.006 


Au 





0.866 ±2E-5 


0.136|3.3E-3 


0.851 ±0.075 




1 


0.151 ±0.004 


1.3E-4|2.7E-3 


0.160 ±0.014 




2 


0.052 ±0.002 


2.1E-3 


0.055 ±0.006 




3 


0.023 ±0.001 


1.7E-3 


0.024 ±0.003 


Ad 





-0.404 ±3E-5 


-0.132|-2.1E-4 


-0.415 ±0.124 




1 


-0.039 ±0.007 


-1.0E-4|-1.7E-4 


-0.050 ±0.022 




2 


-0.011 ±0.002 


-0I-1.3E-4 


-0.015 ±0.009 




3 


-0.004 ±0.001 


Ol-l.lE-4 


-0.006 ±0.005 


Aq 





-0.066 ±0.013 


-0.02|0 


-0.074 ±0.017 




1 


-2.5E-3±1.2E-3 


-3.0E-5 


-0.003 ±0.001 




2 


-3.3E-4±2.0E-4 


0|0 


-4.0E-4±1.0E-4 




3 


-7.0E-5±4.0E-5 


0|0 


-8.0E-5±2.0E-5 


AG 





0.462 ± 0.430 


0.004|1.0E-4 


1.062 ±0.549 




1 


0.079 ± 0.079 


1.0E-5|8.0E-5 


0.184 ±0.103 




2 


0.021 ±0.021 


6.3E-5 


0.050 ±0.028 




3 


0.007 ±0.007 


4.9E-5 


0.017 ±0.010 



Table 4: Moments of the NLO parton densities and their combinations for the present analysis at 
Q^ = 4 GeV^. The value of the respective moment integrating only outside the x-range in which 
currently deep-inelastic scattering data are measured, 0.005 < a; < 0.75, are given for comparison 
(lower|upper part). The errors are the la correlated errors. 
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< X > 


Cp[GeV2] 


C,[GeV2] 


Cp[GeV2] 


Cd[GeV^] 


multiplicative 


additive 


0.060 
0.150 
0.275 
0.425 
0.625 


-0.084 ± 0.245 
-0.229 ± 0.156 
-0.224 ± 0.099 
-0.083 ± 0.140 
0.290 ± 0.417 


0.007 ± 0.287 
0.169 ± 0.431 
-0.226 ± 0.270 
-0.013 ± 0.384 
0.061 ± 1.106 


0.011 ± 0.076 
-0.038 ± 0.037 
-0.045 ± 0.022 
-0.017 ± 0.021 
0.011 ± 0.033 


0.019 ± 0.059 
0.022 ± 0.042 
-0.007 ± 0.025 
0.005 ± 0.027 
0.011 ± 0.031 



Table 5: The higher twist coefficients Cp{x) and Cd{x) as function of x considering both a multi- 
plicative and additive phenomenological ansatz. 
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Figure 1: NLO polarized parton distributions at the input scale Q\ = 4.0 GeV^ (solid line) com- 
pared to results obtained by GRSV (dashed-dotted line) [19], DSSV (long dashed-dotted line) [25j, 
AAC (dashed line) [18], and LSS (long dashed line) [21]. The shaded areas represent the fully 
correlated la error bands calculated by Gaussian error propagation. 
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Figure 2: The polarized parton density xAG{x) at Ql = 4.0 GeV^ as a function of x (solid line). The 
shaded area is the fully correlated Icr statistical error band and the hatched areas are the systematic 
uncertainties. Results from GRSV (dashed-dotted line) [19J, DSSV (long dashed-dotted line) [25], 
AAC (dashed line) [18J, and LSS (long dashed line) [21j are shown for comparison. 
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Figure 3: The polarized parton density xAT.(x) at Ql = 4.0 GeV^ as a function of a; (solid line). The 
shaded area is the fully correlated la statistical error band and the hatched areas are the systematic 
uncertainties. Results from GRSV (dashed-dotted line) [IS], DSSV (long dashed-dotted line) 
AAC (dashed line) [18J, and LSS (long dashed line) [21] are shown for comparison. 
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Figure 4: The spin-dependent structure functions xgl{x), xgf{x) and xg^{x) as a function of x. 
The experimental data are evolved to a common value of Q^ = 5 GeV^. The error bars shown are the 
statistical and systematic ones added in quadrature. The experimental distributions are well described 
(solid curve) within the statistical (shaded areas) and systematic (hatched areas) error bands. The 
curves obtained by GRSV (dashed-dotted) [19J and AAC (dashed) jTH] are shown for comparison. 
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Figure 5: The spin-dependent structure function xgKxjQ"^) as a function of x and Q^. The ex- 
perimental data are compared to the fit result (solid curve) with the statistical error bands (shaded 
areas). The curves obtained by GRSV (dashed-dotted) [19], AAC (short dashed) [H] and LSS (long 
dashed) |21] are shown for comparison. 
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Figure 6: The strong coupling constant a;s(M|) from different DIS measurements, Eqs. (49-56) 
and Refs. [22,25,17]. The yellow band describes the weighted average of a wide range of as{Mz) 
measurements [65j . 
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Figure 7: The additive higher twist coefficients Cp{x), Cd{x) and C„(x) as a function of x. 



28 



References 



[1 

[2 

[s; 

[4] 
[5] 

[6; 

[7 

[9 

[10 
[11 

[12 

[13 
[14 
[15 
[16 
[17 



J. Ashman et al. [European Muon Collaboration], Pliys. Lett. B 206 (1988) 364; Nucl. Pliys. 
B 328 (1989) 1. 

P. L. Anthony et al. [E142 Collaboration], Phys. Rev. D 54 (1996) 6620 



arXiv:hep-ex /96 10007]. 

K. Ackerstaff et al. [HERMES Collaboration], Phys. Lett. B 404 (1997) 383 
|ar Xiv:hep-ex/9703005] . 

K. Abe et al. [E154 Collaboration], Phys. Rev. Lett. 79 (1997) 26 |arXiv:hep-ex/9705012| . 

B. Adeva et al. [Spin Muon Collaboration], Phys. Rev. D 58 (1998) 112001. 

K. Abe et al. [E143 collaboration], Phys. Rev. D 58 (1998) 112003 |arXiv:hep-ph79802357] . 

A. A irapetian et al. [HERMES Collaboration], Phys. Rev. D 75 (2007) 012007 
(a^Xrv:hep-ex/0609039] . 

P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 463 (1999) 339 
|arXiv:hep-ex/9904002| . 

P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 493 (2000) 19 
|arXiv:hep-ph/Q007248| . 

X. Zheng et al.. The JLAB Hall A collaboration, Phys. Rev. C70 (2004) 065207. 

V. Y. Alexakhin et al. [COMPASS Collaboration], Phys. Lett. B 647 (2007) 8 
[arXiv:hep-ex/0609038] . 

K. V. D harmawardane et al. [CLAS Collaboration], Phys. Lett. B 641 (2006) 11 
|arXiv:nucl-ex/0605028] . 

The CLAS collaboration, private communication. 

M. G. Alekseev et al. [The COMPASS Collaboration]. larXiv: 1001.46541 [hep-ex]. 

K. Abe et al. [E154 Collaboration], Phys. Lett. B 405 (1997) 180 |arXiv:hep-ph/9705344| . 



J. Bliimlein and A. Tkabladze, Nucl. Phys. B 553 (1999) 427 [arXiv:hep-ph/9812 478|. 

G. Altarelh, R. D. Ba ll, S. Forte and G. Ridolfi, Acta Phys. Polon. B 29 (1998) 1145 
|arXiv:hep-ph/9803237| ; 

S. Forte, M. L. M angano and G. Ridolfi, Nucl. Phys. B 602 (2001) 585 
arXiv:hep-ph/0101192| . 



[18] Y. Goto et al. [Asym metry Analysis collaboration], Phys. Rev. D 62 (2000) 034017 
|arXiv:hep-ph/0001046]; 

M. Hirai, S. Kumano and N. Saito [Asymmetry Analysis Collaboration], Phys. Rev. D 69 
(2004) 054021 [arXiv:hep-ph/0312112|; 
M. Hirai, S. Kumano and N. Saito, Phys. Rev. D 74 (2006) 014015 |arXiv:hep-ph/0603213l; 



M. Hirai and S. Kumano [Asymmetry Analysis Collaboration], Nucl. Phys. B 813 (2009) 
106 [arXiv:0808.04T3l [hep-ph]]. 

29 



[19] M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, Phys. Rev. D 63 (2001) 094005 
(a^Xw:hep-ph/0011215] . 

[20] C. R. V. B ourrely, J. Soffer and F. Buccella, Eur. Phys. J. C 41 (2005) 327 
|arXiv:hi^^/0502180j . 

[21] E. Leader, A. V. S idorov and D. B. Stamenov, Eur. Phys. J. C 23 (2002) 479 
|arXiv:hep-ph/0111267| . 

[22] E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 80 (2009) 054026 
|arXiv:0908.2390l [hep-ph] ] . 

[23] J. Bliimlein and H. Bottcher, Nucl. Phys. B 636 (2002) 225 ; arXiv:hep-ph/0203155] . 

[24] S. Atashbar Tehran! and A. N. Khorramian, JHEP 0707 (2007) 048 jarXiv: 0705. 26471 [hep- 

ph]]. 

[25] D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. Lett. 101 (2008) 
072001 jarXiv:0804.0422, [hep-ph]]; Phys. Rev. D 80 (2009) 034030 [a rXiv:0904.3821 [hep- 
ph]]. 

[26] B. Adeva et al. [Spin Muon Collaboration], Phys. Rev. D 58 (1998) 112002. 

[27] L. Del Debbio, A. Guffanti and A. Piccione, JHEP 0911 (2009) 060 |arXiv:0907.25 06 [hep- 
ph]]. 

[28] J. Blumlein , H. B5ttcher and A. Guffanti, Nucl. Phys. B 7 74 (2007) 182 
|arXiv:hep-ph /0607200|; Nucl. Phys. Proc. Suppl. 135 (2004) 152 |arXiv:hep-ph/0407089] . 

[29] M. Gliick, E. Reya and C. Schuck, Nucl. Phys. B 754 (2006) 178 [arXiv:hep-ph/0604116]; 
P. Jimenez-Delgado and E. Reya, Phys. Rev. D 79 (2009) 074023 [ arXiv:0810.4274 [hep-ph]]. 

[30] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 64 (2009) 653 
[ar Xiv:0905.3531 , [hep-ph]]. 

[31] S. Alekhin, J. Bliimlein, S. Klein and S. Moch, Phys. Rev. D 81 (2010) 014032 
jarXiv:0908.27M [hep-ph]]. 

[32] A. D. Watson, Z. Phys. C12 (1982) 123; 
W. Vogelsang, Z. Phys. C50 (1991) 275. 



[33] M. Buza, Y. Matiounine, J. Smith et a/., Nucl. Phys. B485 (1997) 420; |hep-ph/9608342] ; 



I. Bierenbaum, J. Bliimlein, S. Klein, (arXiv:0706. 273^1 [hep-ph]] and in preparation. 
[34] S. I. Alekhin, J. Bliimlein, Phys. Lett. B594 (2004) 299. [hep-ph/0404034] . 



[35] M. Lacombe, B. Loiseau, R. Vinh Mau, J. Cote, P. Pires and R. de Tourreil, Phys. Lett. B 
101 (1981) 139; 

W. W. Buck and F. Gross, Phys. Rev. D 20 (1979) 2361; 
M. J. Zuilhof and J. A. Tjon, Phys. Rev. C 22 (1980) 2369; 
R. Machleidt, K. Holinde and C. Elster, Phys. Rept. 149 (1987) 1; 



A. Y. Umnikov, L. P. Kaptari, K. Y. Kazakov and F. C. Khanna, arXiv:hep-ph/9410241 



30 



[36] R. Mertig and W. L. van Neerven, Z. Phys. C 70 (1996) 637 |arXiv:hep- ph/9506451| ; 
W. Vogelsang, Phys. Rev. D 54 (1996) 2023 |a rXiv:hep-ph/9512218] . 

[37] see e.g. : W. Furmanski and R. Petronzio, Z. Phys. C 11 (1982) 293 and references therein; 
G. T. Bodwin and J. W. Qiu, Phys. Rev. D 41 (1990) 2755. 

[38] E. B. Zijlstra and W. L. van Neerven, Nucl. Phys. B 417 (1994) 61, [Erratum-ibid. B 426 
(1994) 245; B 773 (2007) 105]. 

[39] J. Bhimlein and S. Kurth, Phys. Rev. D 60 (1999) 014018 |arXiv:hep-ph79810241|; 



J. A. M. Vermaseren, Int. J. Mod. Phys. A 14 (1999) 2037 arXiv:hep-ph/9806280|. 



[40] J. Bliimlein, Comput. Phys. Commun. 133 (2000) 76 |arXiv:hep-ph/0 003100 1 ; 

J. Bmmlein and S. O. Moch, Phys. Lett. B 614 (2005) 53 |arXiv:hep-ph/0503188 |; 

J. Bliimlein, Comput. Phys. Commun. 180 (2009) 2218 |arXiv:0901.3"T06l [hep-ph]]; 

■arXiv:0901.0837 [math-ph]; 

J. Ablinger, J. Bliimlein, and C. Schneider, in preparation. 

[41] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 48 (1990) 471. 



[42] J. Bliimlein and A. Vogt, Phys. Rev. D 58 (1998) 014020 |arXiv:hep-ph/9712546 |. 

[43] J. Bliimlein and A. Vogt, Phys. Lett. B 370 (1996) 149 |arXiv:hep-ph/9510410 ] ; Phys. Lett. 
B 386 (1996) 350 iarXiv:hep-ph/9606254] . 

[44] J. Bliimlein and A. Vogt, Acta Phys. Polon. B 27 (1996) 1309 [arXiv:hep-ph/960345q] ; 
J. Bartels, B. I. Ermolaev and M. G. Ryskin, Z. Phys. C 72 (1996) 627 
I arXi v:hep-ph/9603204] ; 
Y. Kiyo, J. Kodaira and H. Tochimura, Z. Phys. C 74 (1997) 631 |arXiv:hep-ph/9701365|; 



J. Bliimlein, arXiv:hep-ph/9909449| Lecture Notes in Physics Vol. 546 (Springer, Berlin, 
2000), pp. 42. 

[45] M.J. Alguard et al., E80, Phys. Rev. Lett. 37 (1976) 1261; 
G. Baum et al, E130, Phys. Rev. Lett. 51 (1983) 1135. 



[46 

[47 
[48 
[49 
[50 
[51 
[52 
[53; 
[54; 



L.W. Withlow et al., Phys. Lett B250 (1990) 193. 

M. Arneodo et al., NMC collaboration, Nucl. Phys. B483 (1997) 3. 

K. Abe et al, E143, Phys. Lett. B452 (1999) 194. 

M. Arneodo et al. (NMC), Phys. Lett. B364 (1995) 107. 

D. Adams et al, SMC collaboration, Phys. Lett. B336 (1994) 125. 

K. Abe et al. (E154)), Phys. Lett. B404 (1997) 377. 

K. Abe et al. (E155x)), Phys. Lett. B553 (2003) 18. 

S. Wandzura, F. Wilczek, Phys. Lett. B72 (1977) 195. 

A. Piccione, G. Ridolfi, Nucl. Phys. B513 (1998) 301, |hep-ph/9707478]. 



31 



[55] R. G. Roberts, G. G. Ross, Phys. Lett. B373 (1996) 235, [hep-ph/96012351; 

J. Bliimlein, N. Kochelev, Nucl. Phys. B498 (1997) 285, | hep-ph/9612318 |; Phys. Lett. 
B381 (1996) 296, |h ep-ph/9603397j . 



[56; 

[57 

[58; 
[59; 
[60; 

[61 

[62; 

[63; 

[64 
[65 
[66 

[67 
[68 



[70 

[71 
[72; 



J. Bmmlein and D. Robaschik, Nucl. Phys. B 581 (2000) 449 |arXiv:hep-ph/0002071| . 

J. Bmmlein and D. Robaschik, Phys. Rev. D65 (2002) 096002. |hep-ph/0202 "077l . 

J. Bliimlein, V. Ravindran and W. L. van Neerven, Phys. Rev. D 68 (2003) 114004 
|arXiv:hep-ph/0304292f . 

J. Bliimlein and A. Tkabladze, Nucl. Phys. Proc. Suppl. 79 (1999) 541 
[arXiv:hep-ph/9905524] . 

F. James, CERN Program Library, Long Writeup D506. 

J. Pumplin, D. R. Stump, J. Huston et al, JHEP 0207 (2002) 012, [ hep-ph/020l"T95] . 

A. Airapetian et al. [ HERMES Collaboration ], Phys. Rev. D71 (2005) 012003, 
[hep-ex/0407032; |. 

C. Amsler et al. (Particle Data Group), Phys. Lett. B 667 (2008) 1. 



M. Gliick, E. Reya, A. Vogt, Eur. Phys. J. C5 (1998) 461, |hep-ph/9 806404| . 

S. Bethke, Eur. Phys. J. C 64 (2009) 689 |arXiv:0908.lT35l [hep-ph]]. 

F. D. Aaron et al. [HI Collaboration and ZEUS Collaboration], JHEP 1001 (2010) 109 
[a rXiv:0911.0884 [hep-ex]]. 

S.I. Alekhin, J. Bliimlein, and S.-O. Moch, DESY 10-065, in preparation. 

HI and ZEUS coUab. (V. Radescu et al.). Combined HI and ZEUS Fits Using Low Energy 
Data, talk, DIS 2010, Florence, April 2010. 

R. Abbate, M. Fickinger, A. Hoang, V. Mateu and I. W. Stewart. larXiv: 1004.4894 [hep-ph]: 
larXiv:1006.3080l [hp-ph] . 

T. Gehrmann, M. Jaquier and G. Luisoni, Eur. Phys. J. C 67 (2010) 57 |arXiv:0911.2422l 
[hep-ph]]. 

D. B. Renner. EFXiv: 1002.09251 [hep-lat]. 

H. W. Lin, T. Blum, S. Ohta, S. Sasaki and T. Yamazaki, Phys. Rev. D 78 (2008) 014505 
|arXiv:0802.0"863i [hep-lat]]. 



[73] Y. Aoki et al, Phys. Rev. D 72 (2005) 114505 [arXiv:hep-lat/0411006 |; 
J. Zanotti, private communication for QCDSF; 
T. Yamazaki et a/., Phys. Rev. D 79 (2009) 114505 [arXiv:0904.2039' [hep-lat]]; 

C. Gattringer, C. Hagen, C. B. Lang, M. Limmer, D. Mohler and A. Schafer, Phys. Rev. D 
79 (2009) 054501 [arXiv:0812.1681 [hep-lat]]; 

D. Mohler, private communication for BGR; 
J. Negele, private communication for LHPC; 
T. Korzec, private communication for ETMC. 

32 



[74] Ph. Hagler et al. [LHPC Collaborations], Pliys. Rev. D 77 (2008) 094502 |arXiv:0705.4"295l 
[hep-lat]]; 
J. D. Bratt et al. [LHPC Collaboration] , |ar2£ivilQ0L362Q| [hep-lat] . 

[75] G. Schierholz, private communication. 

[76] M. Gockeler et al., QCDSF collaboration, Phys. Rev. D53 (1996) 2317; Phys. Lett. B414 
(1997) 340; hep-ph/97 11245; Phys. Rev. D63 (2001) 074506; 
S. Capitani et al., Nucl. Phys. (Proc. Suppl.) B79 (1999) 548; 
S. Giisken et al., SESAM collaboration, hep-lat/9901009; 
D. Dolgov et al. [LHPC collaboration and TXL Collaboration], Phys. Rev. D 66 (2002) 



034506 arXiv:hep-lat/0201021 



[77] M. Virchaux and A. Milsztajn, Phys. Lett. B 274 (1992) 221. 

[78] S.I. Alekhin, Phys. Rev. D 68 (2003) 014002 |arXiv:hep-ph/0211096| . 

[79] J. Bliimlein and H. Bottcher, Phys. Lett. B 662 (2008) 336 |arXiv:0802.0408l [hep-ph]]. 

[80] S.I. Alekhin, S. Kulagin and R. Petti, The low-Q deep-inelastic scattering data in the global 
fit of PDFs, Proc. 15th International Workshop on Deep-Inelastic Scattering and Related 
Subjects (DIS2007), Munich, Germany, 2007, Vol. 1, pp. 313. 



33 



